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1 Introduction 



The formulation of the quantitative method of renormalon calculus on the higher level 
of understanding takes its start from the important works of Ref.p^ Ej ; devoted to the 
consideration of e + e~ — > hadrons process, and from the interesting work of Ref. 0, devoted 
to the consideration of deep-inelastic scattering processes. 

After these studies the number of theoretical and practical developments appeared in 
the literature (see reviews of Refs. 

In what is discussed below, we will consider aspects of renormalon calculus, related to 
deep-inelastic scattering (DIS) sum rules. 

It is commonly expected that in the canonical renormalization schemes, say the MS 
scheme, perturbative expansions in small QCD coupling constant a s = a s /(47r) of theo- 
retical expressions for physical quantities, defined in the Euclidean region, are asymptotic 
ones. This means that the difference of the total sums 



D(a) = 1 + E d na n (1) 



n>l 



and their finite sums 



satisfy the following property 



D k (a) = 1 + E d na k (2) 



n=l 



lim a _ r ^0 . (3) 

a k 

In other words the difference between the total series and their finite sum are expressed 
as 

D(a) - D k (a) = 0(a k+1 ) . (4) 

In this case the error of the truncation of the asymptotic series can be estimated by the 
last term of D k (a) , namely d k a k j7j. 

In QCD one expects that in the MS scheme the coefficient function for DIS sum rules, 
normalized to unity can be approximated by the following asymptotic series [B]: 

D k (a) = 1 + £(/3o) k k!(VuV(_i)V + ? (5) 

k>i v ' 

where sign-alternating series with the coefficient i^ uv is generated by the ultraviolet renor- 
malons (UVR), sign-constant asymptotic series with coefficient K 1R result from the con- 
sideration of infrared renormalons (IRR), and a and b are the known numbers, that depend 
from the ratio of the first two coefficients of the QCD /^-function. 

Working within renormalon calculus we will demonstrate that the perturbative and 
non-perturbative contributions to definite DIS sum rules are related. In other words we 
will show that the renormalon approach is working at the boundaries between these two 
regimes in QCD. 

The aim of this article is three fold: 
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• to explain the basic stages of renormalon calculus in QCD, using the simple language; 

• to show that in the asymptotic perturbative expansion of three DIS sum rules, 
namely of the Gross-Llewellyn Smith (GLS), Bjorken- polarized (Bjp) and Bjorken- 
unpolarized (Bjunp) sum rules, may be universal. We will present arguments, based 
on the consideration of the results given in Refs. JD], that these expansions are 
defined by the poles in the closely related Borel images of all three sum rules. 

• We will explain the features, which follow from the consideration of the IRR poles 
in the Borel images of the three DIS sum rules. Moreover our aim is to outline new 
consequences of the IRR calculus. They indicate the existence of relations between 
twist-4 1/Q 2 non-perturbative contributions to the sum rules we are interested in jllj . 
These results form the basis of the new QCD relations between theoretical expres- 
sions for these three sum rules^I], which seem to be supported by the experimental 
data within existing error bars. More critical tests of these relations are proposed. 



2 Renormalon calculus and DIS sum rules 

Let us first express a perturbative QCD series in terms of a Borel integral as 

oo 

D(a 8 ) = J2 d n< (6) 



n=0 

oo 



f°° . . 

= / exp(-5 / (3 a s ) d n —dS 

n=0 n - 

poo 

= / exp(-5/(3 Q a s )B[D](5)d5 , l (7) 
Jo 

where flo = (11/3)Ca — (4/3)7/ iV/ is the first coefficient of the QCD /3-function, with 
Ca = 3, Tf = 1/2, and B[D](S) is the image of the Borel integral. 

At this stage we define the DIS sum rules we will be interested in. The GLS sum rule 
of the uN DIS [12, has the following form 



GLS(Q 2 ) = I [ dx 
2 jo 



Fp(x,Q')+F^(x,Q 



In the Born approximation, this "measures" the number of valence quarks, that are con- 
tained in the nucleon and can thus be considered as the baryon sum rule. In the 
MS scheme, the twist-2 perturbative coefficient function Cgls{Q 2 ) is calculated explic- 
itly, including a 2 s and terms[T3*].|14j. The twist-4 matrix element of the 0(1/Q 2 ) non- 
perturbative contribution to the GLS sum rule is related to the matrix element calculated 
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in Ref.[T5] to be 

where ((0 s )) is defined by the following operator 



uG^ju^u + (u—>d) , 



where 



and 



G 



1 

— ( 

2 



(P\Ol\P) = 2p,((O s )) 



(9) 



(10) 



(12) 



The second sum rules, actively studied both in theory and experiment, is the Bjp sum 
rule^Hl, having the physical meaning of polarized isospin sum rule. Its theoretical 
expression can be defined as 



Bjp(Q 2 ) 



dx 



9a 
6 



Gb ]V {Q 



g l f(x,Q 2 )-g[ n (x,Q 2 
«0 2 » 



Q 2 



o 



1 
Q~ 4 



(13) 



Here Qa = 1-26 is the known /3-decay constant. At the a\ level its perurbative part differs 
from the one of the GLS sum rule by the absence of small "light-by-light" -type terms, 
proportional to the colour structure d abc d abc ^3]. The structure of the power corrections 
to the matrix element of the leading 0(1/Q 2 ) power correction was analytically calculated 
in Ref.|17j. with the useful correction input from the considerations of Ref.|18j. The final 
expressions are presented in a simple-form in the review of Ref.|19j. from which we can 
get: 

1 8 \(U m )) - ^((V m )) , 



«o 2 » 



69 



4 



(14) 



where 



{P,S\U™\P,S) = 2M N S,((U Nii )) 



NS\ 



(P,S\V™\P,S) = 2M N ((V N *)){(S ll P u -S„P,)Ps)} s{ ^ } 



NS\ 



(15) 



and ((U^ s )) and ((^ NS )) are the reduced matrix elements of the local operators from 
Ref . [T7] . namely 



U, 



NS 



V, 



NS 



g s \uG^ v u ~ {u -> d)} 
9s{uG^lsu - {u -> d)} s{V)&} 



(16) 



where S{v, a} stand for symmetrization over the given subscripts and G^ v is defined in 
Eq. (fTTjl . In Ref. [20] the definition of Eq. (14) was used for the estimates of 0(1/Q 2 ) 
corrections to Bjp sum rule, using the three-point function QCD sum rules technique. 
These calculations were then re-analyzed with the same method in Ref. [23] . The numerical 
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results of these calculations will be discussed later. In the work of Ref. [22j a similar analysis 
was done with the help of the same method for the first term in the r.h.s. of Eq. (14), 
while the term, proportional to (M% / 4)((y)) was included into an 0{Mfj/Q 2 ) kinematic 
Al power correction to the Bjp sum rule, which involves the second x 2 moments of the 
leading- twist contribution to g\~ n = g\ — g™ and the twist-3 matrix element, defined 
through the combination of x 2 -weighted moments of the difference of structure functions 
g\~ n and of g p f n = g%(x, Q 2 ) - g%(x, Q 2 ) as 

dt n = £ dxx 2 (2 9 r n (x, Q 2 ) + Sg p 2 ' n (x, Q 2 )) , (17) 

Taking into account this decomposition, it is possible to rewrite a theoretical expression 
for the numerator of the 1/Q 2 contribution, in the way it was done say, in the most recent 
experimental work of Ref. j 



M 2 

^-n = ^ {a V-n + Ad V-n + Af v- n) ? (lg) 

where 

a 2 = f 1 dxx 2 [g{{x,Q 2 )-gl{x : Q 2 )\ (19) 
Jo 

is the target mass correction and 

2m 2 N fr n S, = -4M N S,((U NS )) (20) 
is the twist-4 contribution, which is related to the definition used by us as 

«0 2 » = ^((^ NS )) = -H M Nfi' n ■ (21) 

In other words we have the following relation 

M%fV n = -2«[/ NS » . (22) 

It should be stressed that in the region where the perturbative theory is working well 
enough and the application of the operator-product expansion method is valid (say at 
Q 2 > 2 GeV 2 ), both target mass corrections and twist-3 terms are small and we will 
neglect them in our further considerations 2 . These features were revealed in the process 
of the analysis of Ref. [20 . 

The third sum rule, which was originally derived for purely theoretical purposes, is the 
Bjorken unpolarized sum rule|24j. It can be written down as: 

Fr(x,Q 2 )-Fr(x,Q 2 ) 



Bjunp(Q 2 ) — dx 
Jo 



OWQ 2 )-^-0(i,) • (23) 



2 For completeness we note that there is a minor difference between the the 0(Mjy / 'Q 2 ) coefficients of 
the dxx 2 gi~ n terms in Ref. |2()j and Ref.j2HJ- In the former and latter cases they are equal to (10/9) 
and 1 respectively. 
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It may be also studied in future as the valuable test of QCD both in perturbative and 
non-perturbative sectors. 

As in the previous two cases, the coefficient function CBj U n P (Q 2 ) is calculated up to 
next-to-next-to- leading order a^-corrections|2H],j2ni- The twist-4 matrix element to this 
sum rule was evaluated in Ref.[15j; with the following result: 

«0 3 » = l((O m }) , (24) 
where the matrix element ((0 NS )) is related to the dimension-5 operator 

0^ s = uG^jvjsu - dG^vjsd , (25) 
its matching over nucleon states 

(P\Of\P) = 2^<(0 NS » (26) 

and application of Eq. (j2H) . 

Let us now return to the renormalon calculus. The basic theoretical problem is how 
to define the Borel image B[D}(5) (or the Borel sum) of the integral in Eq. (J7|) for the 
quantities we are interested in. In QCD this problem is usually solved using perturbative 
methods and calculating the corresponding multiloop Feynman diagrams with a one-gluon 
line, dressed by the chains of fermion bubbles (so called renormalon chain insertion). 
These chains are generating sign-alternating asymptotic perturbative series, typical of the 
quantities under consideration, in powers of the expansion parameter Nfa s (where Nf is 
the number of quarks flavours). The contributions of these chains are gauge- invariant, 
but they do not reflect the whole picture of renormalon effects in QCD. The latter begin 
to manifest themselves after application of the naive non-abelianization (NNA) ansatz|2"Tj 
only, namely after the replacement Nf — > — (3/2)/? — Nf — (33/2) in the leading terms of 
the large- Nf expansion. This procedure transforms a large- iVy expansion into a large-/3o 
expansion, which in addition to quark bubbles insertions into the renormalon chain, is 
taking into account the contributions of the gluon- and ghost-bubbles insertions as well 
(though neglecting definite one-loop insertions into the gluon-quark-antiquark vertex, 
which should be also considered in the process of rigorous calculation of the coefficient 
(3q). The application of the NNA approach allowed the authors of Ref. |2E] to formulate 
the extension to higher orders of the BLM-approach [29J. Technically, the work of Ref. 
[2~%] supports the results of the first successful formulation of the BLM procedure to the 
next-to-next-to-leading order [3U] Moreover, these two works pushed ahead the study of 
the BLM procedure in higher orders [31] . In principle, the relations of the results of Refs. 
[3TH |2E1 EH] need more detailed considerations. In view of the lack of space we will avoid 
discussions of this subjects here. 

The Borel images calculated by this procedure for the GLS and Bjp sum rules coincide 
and have the following form [5] : 

B[C BiP ] W = B[Coi,s]ffl = - ■ < 27 > 
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They contain the IRR poles at 5 = 1 and 5 = 2 and the UVR poles at 5 = — 1 and 5 = —2. 
Note that the 5 = — 1 UVR poles in Eq. (J27j) are suppressed by a factor (l/2)exp(— 10/3) = 
0.018, relative to the dominant IRR poles at 5 = 1 10 . Therefore, in the asymptotic 
structure of the perturbative QCD effects in the expressions for Cgls(Q 2 ) ~ CBj P (Q 2 ) 
(where we neglect the small "light-by-light-type" effects, contributing to Cgls(Q 2 )) the 
sign-constant part in Eq. (jSJ) dominates strongly with respect to the sign-alternating 
contribution, generated by 5 = — 1 UVR. The scheme-dependence of these results are not 
so obvious, Indeed, the suppressions of 5 = 1 UVR with respect to 5 = 1 IRR is related to 
the application of the MS-scheme which we are using throughout the whole work. In fact 
in this scheme the IRR renormalons are not suppressed. However, there is the procedure, 
when the situation is reversed- the IRR are absent, but UVR may exist. This feature is 
manifesting itself for the models with "frozen" coupling constant (see e.g. [33] )• 

Returning to the large- Nf expansion of the perturbative expressions 

f~i oo 

C Bjp (Q 2 ) = C GLS (Q 2 ) = l + -^-E r n (T f N f a s r , (28) 



TfNf 



where C F = 4/3, T f = 1/2 and 



4 d 



-l 



r n = lim^---j B[C m ](5) , (29) 

we arrive at the following expansion in powers of x = TfNfa s , namely 

2 920 , 38720 4 . , 

J2r n x n = -3x + 8x 2 -— x 3 + -^-x 4 + ... , (30) 

which is known in the MS scheme up to 0(a 9 s Nj) terms(5]. Using now the traditional 
MS-scheme expansion in terms of the orders in a s /ir = 4a s , one can compare the results 
of explicit perturbative calculations of 

^B jP (Q 2 ) = l + E^(v) n ( 31 ) 

with the known numbers 

n = -l (32) 

r 2 = -4.5833 + 0.33333V/ (33) 
r 3 = -41.440 + 7.6073JV/-0.17747JV? (34) 

obtained at 0(a 2 s ) in Ref-fTS] and at 0(a^) in Ref. [Tl]. with the results of the application 
of the NNA procedure[2ZI to the estimates of the perturbative QCD corrections from large- 
Nf expansion of Eq. (J30|) 3 . Performing the shift Nf — > Nf — 33/2 in the second, third 
and fourth terms in Eq. ()30|). we arrive at the following estimates in the MS scheme [Til] : 

r NNA = -5.5 + 0.33333^ (35) 

r NNA = _48.316 + 5.8565V/ - 0.17747A^| (36) 

r NNA = _466.00 + 84.728V/- 5.1350 ^ + 0.10374A^ 3 . (37) 



3 It is worth noting that similar NNA analysis was performed previously, in Ref . |32| . for the e + e 
annihilation Adler D-function. 
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Reasonable agreement can be observed between the sign structure and values of the NNA 
estimates and the results of explicit calculations (compare the estimates of Eqs. (|35|) and 
(jnHJ) with the numbers in Eqs. (|3*3*|) and (jSjj) , respectively). As to the prediction for r^ NA , 
it may serve as a guide for understanding the rate of growth of the coefficients of the 
perturbative series generated by the single renormalon-chain approximation. 

Consider now the Bjunp sum rule, which is defined in Eq. (|23)1 . Within the large- Nf, 
expansion its perturbative coefficient function 



n>l 



^B j punp(g 2 ) = i + E^( v vJ (38) 



was calculated in the MS scheme and large- Nf expansion up to a Q(q.A^|)-terms[TU] . 
Following the logic of our work, we present here the results for the first 4 terms only: 

64 2 2480 , 113920 4 , , 

Y f nX n = -2x + —8x 2 x 3 + — x 4 + . . . 39 

V 9 81 729 v ; 

As was already mentioned above, the explicit results of calculations of the perturbative 
contributions to the Bjunp sum rule 

C Bj unp(Q 2 ) = l+E^(-)" (40) 
n>l 

are known up to the order 0(a 3 ) level. These results are: 

h = -2/3 (41) 
f 2 = -3.8333 + 0.29630^/ (42) 
f 3 = -36.155 + 6.3313JV)- 0.15947 Nf (43) 



where f 2 was calculated in Ref.[25J while f 3 was evaluated in Ref. |26j. Applying now the 
NNA procedure to the series of Eq. (|39|h we find that, in the MS scheme, the estimated 
coefficients of the Bjunp sum rules have the following form 10J: 

~nna = _4.8889 + 0.29630^/ (44) 
~nna = _43_ 414 + 5.2623^/ - 0.15947^ (45) 
~nna = _ A57 Q2 + 83.094^/ - 5.0360 JVj + O.lOUANf . (46) 



The estimate of Eq. (|44j) is in agreement with its exact partner of Eq. (|42|). The same 
situation holds for the 0(a 3 ) corrections (compare Eq. (|43|) with Eq. It should be 

stressed, that the similarity of the next-to-next-to- leading-order MS-scheme perturbative 
QCD contributions to the Bjp and Bjunp sum rules was previously noticed in Ref . [3"4~j . 
although no explanation of this observation was given. Now, within the NNA procedure, it 
is possible to generalize this observation to higher-order level. Indeed, the NNA estimates 
of the 0(a 4 ) corrections to the Bjp and Bjunp sum rules have a similar expressions as well. 
These facts may indicate the close similarity in the full perturbative structure 
of the QCD corrections to the Bjunp sum rule, the Bjp sum rule and the GLS 
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sum rule (provided the "light-by-light-type" terms will not drastically modify the values 
of perturbative terms in the latter case in the one-renormalon chain approximation). Note 
that, generally speaking, from this order of perturbation theory the diagrams from the 
second renormalon chain are starting to contribute to the quantities under consideration. 
These diagrams may influence the asymptotic behavior of the the series considered j35|. In 
view of this it seems that it is more rigorous to use, in the phenomenological application, 
the order of a^-terms, estimated in Ref.[36 3 using the PMS approach^?] and the effective- 
charges approach, developed in Ref.j3H]. However, since in this work we concentrated 
ourselves on the structure of the QCD expressions, obtained in the one-renormalon chain 
approximation, we will avoid more detailed discussions of the possible influence of the 
multi-renormalon chain contributions to the results of our studies. 

The observed in Ref. (3T similarity of the next-to-next-to-leading-order approximations 
for the Bjp and Bjunp sum rules was attributed in Ref.fTU] to the fact that the dominant 
6 = 1 IRR contribution to the Borel images of these sum rules enters with identical 
residues. Indeed, the Borel images in the Borel integrals of Eq. (7) for the Bjunp and Bjp 
sum rules turn out to be closely related [TU]. namely 

^■(^)^- tl -T^) ' ■ (47) 

Comparing Eq. (27) with Eq. (jTfjl one can convince oneself that the residues of the 
poles at 5 = 1 in these two expressions are really the same and are equal to the factor 
-(8/3)exp(5/3). 

Notice also the absence of 5 = —1 UVR pole and the existence in Eq. (pTTj) of a 
S = — 2 UVR pole together with the leading 5=1 IRR one. Thus we are observing 
one more interesting fact: the structure of the Borel image for the Borel sum, related 
to the Bjunp sum rule, is dual to the structure of leading renormalon contributions to 
the Borel image of the Borel sum for the e + e~ annihilation Adler D-function. Indeed, 
in the latter case the leading IRR is manifesting itself at 5 = 2, while the leading UVR 
pole is appearing at 5 = — 1 (the general structure of renormalon singularities in the e + e~ 
annihilation channel was analyzed in Ref. pQ, while the concrete MS-scheme calculations 
of the corresponding Borel image were done later on in Refs.[39j and |40j). 

The absence of 5 = 1 IRR in the Borel sum of the e + e~ annihilation channel is related 
to the absence of 0(A 2 /Q 2 ) non-perturbative power correction in the standard variant of 
the operator product expansion formalism, applied to the theoretical expression for the 
e + e~ annihilation Adler D-function. Indeed, the existence of lowest dimension-4 quark and 
gluon condensates |H] in this channel can be associated in terms of renormalon language 
with the existence of the leading IRR pole, which in case of "Borelization" of the Adler D- 
function is appearing at 5 = 2. However, as was already discussed above, the dimension-2 
non-perturbative corrections enter into the theoretical expressions for the three DIS sum 
rules we are interested in. In the IRR language, this corresponds to the appearance of 
a 5 = 1 IRR polejH], which manifests itself in the concrete results of Refs.p].[TU] (see 
Eqs. (27) and (47)). Thus, it should be stressed that the structure of singularities of the 
Borel sums (or images) is not universal and depends from the physical quantity under 
consideration. 



3 IRR for DIS sum rules and the values of twist-4 
corrections 



In addition to controlling the sign-positive n\ growth of the asymptotic series the existence 
of 5 = 1 IRR gives an ambiguity in taking the Borel integral of Eq. (7) over this pole. 
In the case of large (3q expansion and for the series we are interested in, this ambiguity 
was estimated in Ref. [S]. Moreover, 5 = 1 IRR generates the negative power suppressed 
correction which has the following expression: 



AC. 



sum rules 



32exp(5/3) A- 



MS 



9/3o Q 2 



(48 



Notice, that it has the same negative sign as the residue of 5 = 1 IRR. 



This estimate may be coordinated with the definition of the twist-4 matrix element 
in the sum rules we are interested in. Therefore, we will make the assumption that 
the identical values and signs of the IRR induced power-suppressed term indicate that 
the values of twist-4 contributions to the expressions of GLS, Bip and Bjunp sum rules, 
normalized to unity, should have the same negative sign and a similar closed value [TT] . 
This assumption is similar to the known in the renormalon-oriented literature guess of 
"universality" 



Let us check this assumption, considering the following expressions for the sum rules 
we are interested in 



GLS(Q 2 ) 
Bjp(Q 2 ) 
Bjunp(Q 2 ) 



A 



9A 

6 



1 - 4a. - 0(a 2 s ) - ^ 
1 - 4a s - 0(a 2 



Q 2 



B 

Q 2 



1 - -a a - 0(a 2 ) - 



_C_ 

Q 2 



(49) 
(50) 
(51) 



where A = ((Oi))/3, B = ((O2)) (6/(74) and C = ((O3)) and compare in Table 1 the results 
of different theoretical and phenomenologically based evaluations of the twist-4 parameters 
A, B and C. 

In the case of the GLS and Bjunp sum rules the results of the original application of 
the three-point function QCD sum rules method gave ((0 s )) = 0.33 GeV 2 and ((0 NS )) = 
0.15 GeV 2 , with over 50% error barsjlH], while the three-point function estimates for the 
modified results of calculations of the twist-4 parameter of the Bjp sum rule resulted in 
the following value M^/f - ™ = —0.18 ± 0.09 GeV 2 in the region where nucleon target 
mass corrections of 0(Mff/Q 2 ) and twist-3 contribution may be neglected [2U]. As was 
already mentioned, these calculations were re-examined using three-point function QCD 
sum-rules approach in Refs. [2*T] and [22] • In the first case the obtained result turned 
out to be larger than the original results from Ref.[2*U] and has the following value of 
M^/f"™ = —0.634 ±0.317 GeV 2 [21], while in the latter case it was considerably smaller, 
namely M|r/f~ n = —0.07 ± 0.035 |22j, although within 50% theoretical uncertainty we 
adopt for all calculations within three-point function QCD sum-rules approach, this value 
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Table 1: The results for twist-4 contributions to the GLS, Bjp and Bjunp sum-rule ex- 
pressions of Eqs. (49)-(51). 





A [GeV 2 ] 


B [GeV^] 


C [GeV^] 


QCD sum rules (Ref.03|) 
QCD sum rules (Ref.fSg 
QCD sum rules (Ref.j2H 
QCD sum rules (Ref. [22]) 
Instanton model (Ref.[46j) 
Instanton model (Ref.jUj) 
Experiment (Ref.|48j) 
Experiment (Ref. j50j) 


0.098 ± 0.049 
0.158 ±0.078 
0.078 ±0.039 

0.04 ±0.13 


0.063 ±0.031 
0.223 ±0.118 
0.025 ±0.012 
0.087 ±0.043 

0.098 ±0.028 


0.133 ±0.065 
0.16 ±0.08 

0.16 ±0.08 



does not disagree with the results obtained in Ref. |20j. The relatively high difference 
between the central values of estimates of Refs. [21] and [201 is explained by the fact that in 
the former analysis the additional corrections to the perturbative side of the corresponding 
QCD sum rules are included and the continuum term to the nucleon pole of the low-energy 
side of this sum rule is explicitly retained. This leads to better stability of the extracted 
value of matrix element with respect to the Borel parameter and increases its central value. 
Note, however, that the theoretical error of the three-point function of the QCD sum rules 
result of Ref. [21] is considerably underestimated. We fix it as 50 % uncertainty, which to 
our point of view is typical to all three-point function QCD sum rules results. 

In Table 1 we present the estimates of twist-4 corrections to different DIS sum rules, 
obtained with the help of the three-point function QCD sum-rules approach and compare 
them with the results of the application of different theoretical approach, based on the 
picture of the QCD vacuum as a "medium" of instantons[44 . This picture was further 
developed in the method in Ref. [13] and applied for estimating twist-4 contributions to the 
GLS sum rule and Bjp sum rule in Ref. |46| . while the number for the twist-4 contribution 
to the Bjunp sum rule, which follows from this approach, was presented in Ref. |47j . In 
the absence of estimates of theoretical uncertainties within this approach, we will apply 
to them the careful 50% estimate as well. All these results support the original results 
of the three-point function QCD sum rules calculations of the twist-4 corrections to the 
GLS, Bjunp sum rules|I3] and Bjp sum rulej2D], though the additional three-point function 
QCD sum rules cross-check of the results of Ref. [20] m ay be rather useful. 

The experimentally motivated value of the twist-4 contribution to the Bjorken sum 
rule M^/f~ n = —0.28 ± 0.08 GeV 2 [IH] was obtained by means of integrating in x the 
numerator of the dimensionless h(x)/Q 2 contributions, extracted from the fits of world 
average data for g{(x, Q 2 ) and g™(x, Q 2 ) performed in Ref. [1H1- From the results of Table 
1 one can see that the agreement with the QCD sum-rules calculations of Ref. [20] and 
instanton-based calculations of Ref. [1H] is more than qualitative. 

The experimentally inspired estimate for the twist-4 contribution to the GLS sum rule 
was obtained only recently [SU] as a result of the integration of ^-dependence of the twist-4 
contribution h(x) / Q 2 , extracted in the works of Ref. [ST] devoted to the analysis of xF 3 data 
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of CCFR collaboration. One can see that the central value of the contribution is negative 
(in fact it comes with the negative sign in the sum rule) but has rather large uncertainties. 
So, at the present level we cannot obtain from this estimate even qualitative information 
and additional work on its improvement is needed. 

To conclude, we present the final results for the GLS, Bjp and Bjunp sum rules, where 
for definiteness the twist-4 matrix elements are estimated using the central values of the 
three-point function QCD sum-rules results from Refs. |4H] and [2U] : 



GLS(Q 2 ) 
Bjp(Q 2 ) 
Bjunp(Q 2 ) 



9a 
6 



1 -4a s -0(a 2 s ) 
1 - 4a s - 0(a 2 s 



0(a 2 s ) 



0.098 GeV 2 
Q 2 

0.063 GeV 2 



Q 2 

0.133 GeV 2 



Q 2 



(52) 
(53) 
(54) 



It should be stressed that they all have the same negative sign and within existing 
theoretical uncertainties are in agreement with each other. This fact was anticipated by 
the identical value of the ambiguity, generated by the 6 = 1 IRR pole of the Borel images of 
all these three sum rules (see Eq. (48)). Moreover, as follows from the results of application 
of the single-renormalon chain approximation in the perturbative sector presented in Sec. 2, 
we may expect a similar asymptotic behavior of the perturbative corrections to all these 
three sum rules (compare Eqs.(33)-(36) with Eqs.(42)-(45)). It is interesting that the 
similar property is manifesting itself in perturbative series under investigations at the 
O(o^) level, studied within scheme- invariant approaches in Ref.[36j. 

These facts give us the idea that the sum rules we are interested in are closely related 
and that, in the region, where we can neglect target mass corrections and twist-3 con- 
tributions to the Bjp sum rule and quark-mass dependent corrections (say in the region 
Q 2 > 2 GeV 2 ) we can write down the following basic relation [TT] : 

Bjp(Q 2 ) « (W18)GLS(Q 2 ) « (W6)Bjunp(Q 2 ) . (55) 



In the next section we will present more detailed considerations of the experimental 
consequences of these relations then those, that are briefly outlined in Ref.|llj. 



4 IRR- inspired relations and experiment 

In order to test whether our basic relation Eq. (J55|) is respected by experiment, we first 
present the results of the extraction of the GLS sum rule by combining CCFR neutrino 
DIS data with the data for other neutrino DIS experiments for 1 GeV 2 < Q 2 < 15 GeV 2 
|52j . It is known that the weighted extraction of a s (Mz) from these data result in the 
rather rough value a s (Mz) = 0.115±o < j ) 2 9 , which, is in agreement with a s (Mz) = 0.115 ± 
0.001 (stat) ±0.005 (syst) ±0.003 (twist) ± 0.0005 (scheme), extracted in Ref.jSSl from 
the previous CCFR data for the GLS sum rule at Q 2 = 3 GeV 2 |Hlj. However, for our 
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Table 2: The results for the GLS sum rule from Ref. \5H$ 



Q 2 [GeV 2 ] 


GLS sum rule 


2.00 


2.49 ±0.08 ±0.14 


3.16 


2.55 ±0.08 ±0.10 


5.01 


2.78 ±0.06 ±0.19 


7.94 


2.82 ±0.07 ±0.19 


12.59 


2.80 ±0.13 ±0.18 



Table 3: The comparison of the results of application of Eq. 155]) with direct experimentally 
motivated numbers 



Q 2 [GeV 2 ] 


Bjp from Table 1 


Bjp SR (exp) 


2.00 


0.174 ±0.006 ±0.010 


0.169 ±0.025 [Ref. 


3.16 


0.178 ±0.004 ±0.007 


0.164 ±0.023 [Ref. 


5.01 


0.195 ±0.004 ±0.013 


0.181 ± 0012 (stat) ± 0.018 (syst) [Ref.|5E|] 


7.94 


0.197 ±0.005 ±0.013 




12.5 


0.196 ±0.009 ±0.013 


0.195 ±0.029 Ref. jHZj 



purposes we will not need to re-extract a s (M z ) values from the GLS sum rule results of 
Ref. |52j . but will use these, which are presented in Table 2. 

To estimate the values of the Bjp sum rule from the results of Table 2 we will use our 
main equation (]53]l and will compare them with available experimental data for the Bjp 
sum rule. The results of these studies are presented in Table 3. 

One can see that though the central values of estimated numbers for the Bjp SR are 
higher than the results of the SLAC E143 collaboration they agree within error bars. 
It is also interesting to compare the result from Table 3 with the value of the Bjp sum 
rule extracted in Ref.jSHl from the SLAC and SMC data Bjp(3 GeV 2 ) = 0.177 ±0.018 and 
which, within error bars, do not contradict the value Bjp(3 GeV 2 ) = 0.164 ± 0.011 used 
in the work of Ref. [59 . It is rather inspiring that within error bars these results agree 
with the GLS sum rule value at Q 2 = 3.16 GeV 2 . The same feature holds for the Bjp sum 
rules at Q 2 = 5 GeV 2 , namely for the SMC result of Ref. |56j . Thus we think that within 
existing uncertainties our approximate IRR-inspired basic equation (jST)]) is supported by 
existing experimental data. 

5 Conclusions 

We demonstrated that the existing phenomenological data do not contradict the basic re- 
lation of Eq. (|55j) and therefore the reliability of the one-renormalon chain approx- 
imation of the theoretical quantities under consideration. For its more detailed 
studies, we may rely on the appearance of Neutrino Factory data for all sum rules, which 
enter in Eq. (|55p. In fact it may provide rather useful data not only for the GLS and 
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Bjp sum rules, but for the Bjunp sum rule as well (for a discussion of this possibility see 
Refs.[Ul,|I>I]). 

Another interesting option of the relation of Eq. (|55|) is to analyze the sources of its 
possible violation in the lower energy region of over Q 2 m 1 GeV 2 , where one may compare 
the CCFR data for the GLS sum rule at the energy point Q 2 = 1.26 GeV 2 |S21 and the 
JLAB data for the Bjp sum rule at Q 2 = 1.10 GeV 2 P3j. 

To conclude this section, we would like to emphasize that the problems considered 
by us in this work are complementary to the considerations of Ref.[62j. In the former 
analysis, the GLS and Bjp sum rules were determined in the high energy point of over 
Q 2 = 12.33 GeV 2 from the generalized Crewther relation constructed in using 
the extension of the BLM approach of Ref. and the analysis of e + e~ annihilation data 
from Ref.jOl]. Certainly, the renormalon- chain insertions are absorbed in this approach 
into the BLM scale. However, the considerations within this language of the high-twist 
effects is still missed. It may be of interest to think of the possibility of evaluating high- 
twist contributions to the Crewther relation, which relates, in the Eucledian region we 
are working massless QCD perturbative contributions to the Adler D-function of e + e~- 
annihilation with the perturbative corrections to the GLS and Bjp sum rules. 
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